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1.0  INTRODUCTION 

The  present  tests  were  conducted  to  obtain  heat-transfer-rate  data  on  the  Space  Shuttle 
Integrated  Vehicle  at  simulated  ascent  flight  conditions.  The  components  of  the  vehicle 
produce  a  complex  flow  field  with  many  shock  interactions.  Although  progress  has  been 
made  in  analyzing  shock  interference  heating  (Ref.  1),  experimental  data  are  still  relied 
upon  to  provide  vehicle  design  information.  One  technique  resulting  from  analytical  studies 
is  the  ratioing  of  interference  heating  to  undistrubed  values.  During  the  present  tests,  the 
external  tank  was  tested  alone  so  that  these  heating  ratios  could  be  constructed  and  also 
to  confirm  the  data  quality. 

The  tests  were  conducted  in  the  Supersonic  Wind  Tunnel  (A)  of  the  von  Kkrmkn 
Gas  Dynamics  Facility  (VKF)  at  Mach  numbers  2.5,  3.5,  4.5,  and  5.5  and  a  free-stream 
unit  Reynolds  number  of  5  million  per  foot.  Angle  of  attack  was  varied  from  -5  to  10 
deg  and  sideslip  angles  varied  from  0  to  -6  deg. 

2.0  APPARATUS 


2.1  WIND  TUNNEL 

Tunnel  A  is  a  continuous,  closed-circuit,  variable  density  wind  tunnel  with  an 
automatically  driven  flexible-plate-type  nozzle  and  a  40-  by  40-in.  test  section.  The  tunnel 
can  be  operated  at  Mach  numbers  from  1.5  to  6  at  maximum  stagnation  pressures  from 
29  to  200  psia,  respectively,  and  stagnation  temperatures  up  to  750'‘R  (M^  =  6).  Minimum 
operating  pressures  range  from  about  one-tenth  to  one-twentieth  of  the  maximum  at  each 
number.  The  tunnel  is  equipped  with  a  model  injection  system  which  allows  removal  of 
the  model  from  the  test  section  while  the  tunnel  remains  in  operation.  A  schematic  view 
of  Tunnel  A  is  presented  in  Fig.  1. 

2.2  MODEL 

The  test  model  was  a  0.0175-scale,  thin  skin,  thermocouple-equipped  model  of  the 
Rockwell  International  Vehicle  5  configuration  of  the  Space  Shuttle.  Rockwell 
International  fabricated  the  model  and  supplied  the  model  drawings.  A  sketch  of  the  model 
showing  the  model  coordinate  system  and  the  overall  full-scale  dimensions  is  presented 
in  Figure  2.  The  integrated  model  was  composed  of  the  orbiter  vehicle,  external  tank, 
and  two  solid  rocket  booster  motors  that  are  identified  in  Fig.  3. 

The  external  tank  was  designed  with  two  interchangeable  nose  tip  configurations. 
The  two  nose  tip  profiles  are  shown  in  Fig.  4.  The  smooth  nose  configuration  in  Fig. 
4a  is  an  extension  of  the  ogive  contour  of  the  forward  section  of  the  external  tank.  The 
nipple  nose  configuration.  Fig.  4b,  has  a  spherical  nose  tip  extended  from  the  ogive  contour. 
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The  model  was  constructed  of  17-4  PH  stainless  steel.  The  nominal  skin  thickness 
was  0.030  in.  at  the  instrumented  areas.  All  thermocouples  were  spot-welded  to  the  inner 
surface  of  the  model  and  the  skin  thickness  was  measured  at  each  thermocouple  location. 

Boundary-layer  trips  were  designed  to  be  form-fitted  to  the  orbiter  and  the  solid 
rocket  boosters.  The  orbiter  and  solid  rocket  booster  trips  consisted  of  0.020-in.-diam 
stainless  steel  balls  spot-welded  to  a  thin,  stainless  steel  band.  The  trips  attached  to  the 
orbiter  at  an  Xq/Lq  of  0.04  and  on  the  solid  rocket  boosters  at  an  Xs/L*  of  0.019.  The 
trip  ring  on  the  external  tank  consisted  of  a  band  of  No.  36  Carborundum®  grit  located 
at  an  xt/Lt  of  0.069.  A  pictorial  view  of  the  trip  ring  on  the  external  tank  with  the 
smooth  nose  attached  is  shown  in  Fig.  5. 

The  model  configurations  tested  consisted  of  the  external  tank  alone  with  either  the 
smooth  nose.  Fig.  5,  or  the  nipple  nose.  Fig.  6,  installed  and  all  protuberances  removed 
and  the  integrated  model  with  the  nipple  nose  installed  on  the  external  tank.  Fig.  3. 
No  data  were  recorded  for  the  orbiter  model  during  this  test. 

2.3  INSTRUMENTATION  AND  MEASUREMENT  PRECISION 

Tunnel  A  stilling  chamber  pressure  is  measured  with  a  15-,  60-,  150-,  or  300-psid 
transducer  referenced  to  a  near  vacuum.  Based  on  periodic  comparisons  with  secondary 
standards,  the  uncertainty  (a  bandwidth  which  includes  95  percent  of  the  residuals)  of 
these  transducers  is  estimated  to  be  within  ±0.2  percent  of  reading  or  ±0.015  psia, 
whichever  is  greater.  Stilling  chamber  temperature  is  measured  with  a  copper-constantan 
thermocouple  with  an  uncertainty  of  ±3“F  based  on  repeat  calibrations. 

The  model  temperatures  were  measured  with  iron-constantan  thermocouples  with  an 
estimated  uncertainty  of  ±0.5  percent.  There  were  837  thermocouples  installed  on  the 
model.  A  Beckman®  210  analog-to-digital  converter  was  used  in  conjunction  with  a  CDC 
160-A  computer  to  record  the  temperature  data. 

The  standard  Tunnel  A  instrumentation  system  is  designed  to  record  a  maximum 
of  97  thermocouples.  To  more  efficiently  accommodate  the  837  thermocouples  on  this 
test  required  the  installation  of  additional  equipment.  A  new  terminal  board  was  fabricated 
that  allowed  291  thermocouples  to  be  connected  at  one  time.  An  existing  three-position 
switch  was  incorporated  that  was  used  to  select  a  group  of  97  thermocouples  to  be  recorded 
during  each  model  injection.  All  the  model  thermocouple  leads  were  terminated  at 
quick-disconnect  plugs.  These  plugs  were  easily  interchanged  with  mating  plugs  that 
attached  to  the  terminal  board.  In  this  manner,  a  preselected  group  of  291  thermocouples 
could  be  rapidly  connected,  and  by  varying  this  hookup  it  was  possible  to  record  data 
from  all  837  thermocouples. 
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3.0  PROCEDURE 

3.1  TEST  CONDITIONS 

The  tests  were  conducted  at  the  following  nominal  conditions: 


hREF  . 


M 

oo 

Po,  psia 

To.'R 

Btu/sec-ft^-’R 

Re  ,  ft-’ 

oo' 

2.5 

32 

660 

0.062 

5.0  X  106 

3.5 

57 

680 

0.055 

5.0  X  106 

4.5 

95 

665 

0.047 

5.0  X  106 

5.5 

178 

740 

0.045 

5.0  X  106 

Data  were  obtained  at  angles  of  attack  of  0,  ±5,  and  10  deg  and  yaw  angles  of  0,  3, 
and  6  deg.  A  complete  test  matrix  with  the  model  configurations  tested  is  included  in 
Table  1. 

Uncertainties  of  the  basic  tunnel  parameters  were  estimated  from  repeat  calibrations 
of  the  po  and  To  instruments  and  from  the  repeatability  and  uniformity  of  the  tunnel 
flow  during  calibrations.  The  parameters,  Po,  To,  and  with  their  uncertainties,  were 
then  used  to  compute  the  uncertainties  in  the  other  parameters  dependent  on  these  by 
means  of  the  Taylor  series  method  of  error  propogation. 

Approximate  uncertainties  in  tunnel  flow  parameters  are  shown  as  follows: 


Uncertainty  f±).  percent 


M 

OB 

M 

OO 

Po 

To 

Re 

oo 

2.5 

0.8 

0.2 

0.5 

1.3 

3.5 

0.5 

0.2 

0.5 

1.3 

4.5 

0.4 

0.2 

0,5 

1.3 

5.5 

0.3 

0.2 

6.5 

1.3 

3.2  TEST  PROCEDURE 

The  desired  model  attitude  was  established  prior  to  injecting  the  model  into  the  flow. 
When  the  model  reached  tunnel  centerline,  the  model  was  immediately  translated  forward 
into  the  undistrubed  tunnel  flow.  The  thermocouple  outputs  were  scanned  approximately 
20  times  per  second  starting  prior  to  model  injection  into  the  airstream  and  continuing 
about  4  sec  after  the  model  reached  the  tunnel  centerline.  After  each  injection,  the  model 
was  cooled  to  an  isothermal  state  using  high-pressure  air. 


7 


AEDC-TR-7S-153 


3.3  DATA  REDUCTION  AND  PRECISION 

I 

The  reduction  of  thin  skin  thermocouple  data  normally  involves  only  the  calorimetric 
heat  balance  which  in  coefficient  form  is: 

dTw/dt 

h(To)  -  wbcp  (1) 

Radiation  and  conduction  losses  are  neglected  in  this  heat  balance,  and  data  reduction 
simply  requires  evaluation  of  dTw/dt  from  the  temperature-time  data  and  determination 
of  model  material  properties.  For  the  present  tests,  radiation  effects  were  negligible; 
however,  conduction  effects  can  be  significant  in  several  regions  of  the  models.  To  permit 
identification  of  these  regions,  and  to  improve  evaluation  of  the  data,  the  following 
procedure  was  used. 


Separation  of  variables  and  integration  of  Eq.  (1)  assuming  constant  w,  b,  Cp,  and 
To  yields 


h(To) 
w  b  Cp 


Differentiation  of  Eq.  (2)  with  respect  to  time  gives 

MTq)  d  r  To  -  Twq 

w  b  Cp  “  "dr  -Tw  J 


(2) 


(3) 


Since  the  left  side  of  Eq.  (3)  is  a  constant,  plotting  versus  time  will  give  a  straight 
line  if  conduction  is  negligible.  Thus,  deviation  from  a  straight  line  can  be  interpreted 
as  a  conduction  effect. 


The  data  were  evaluated  in  this  manner,  and  generally  a  linear  portion  of  the  curve 
was  used  for  all  thermocouples.  A  linear  least-squares  curve  fit  of  fin[(To  -  Twj)/(To  - 
Tw)]  versus  time  was  applied  to  the  data.  The  data  from  a  typical  thermocouple  are 
presented  in  Fig.  7.  The  value  of  the  initial  slope  is  high  when  the  model  is  located  on 
the  tunnel  centerline  in  the  rear  test  section.  As  the  model  is  translated  forward  out  of 
this  region  into  the  undisturbed  tunnel  flow  a  new  linear  slope  is  established  that  is  the 
desired  value  of  d/dt  fin[(To  -  Twi)/(To  -  Tw)].  Generally,  for  all  thermocouples  with 
a  value  of  xt/Lt  less  than  or  equal  to  0.05  the  curve  fit  was  started  at  the  time  the 
model  reached  the  tunnel  centerline.  The  curve  fit  for  thermocouples  with  a  value  of 
Xt/Lt  greater  than  0.05  and  less  than  or  equal  to  0.1  started  1.0  sec  after  the  model 
reached  centerline.  The  curve  fit  for  the  remaining  thermocouples  started  2.0  sec  after 
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centerline  to  ensure  that  the  model  was  clear  of  any  disturbance  in  the  aft  part  of  the 
test  section.  The  curve  fit  extended  for  a  time  span,  which  was  a  function  of  the  heating 
rate,  as  shown  below: 


_ Range _  Number  of  Points 

dTw/dt  >  32 
16  <  dTw/dt  ^  32 

8  <  dTw/dt  ^  16 

4  <  dTw/dt  ^  8 

2  <  dTw/dt  ^  4 

1  <  dTw/dt  ^  2 

dTw/dt  <  1 

The  above  time  spans  were  adequate  to  keep  the  evaluation  of  the  right  side  of 
Eq.  (3)  within  the  linear  region.  The  linearity  of  the  fit  was  substantiated  by  visual 
inspection  of  the  cases  in  question.  Strictly  speaking,  the  value  of  Cp  for  the  material 
was  not  constant,  as  assumed,  and  the  following  relation: 

Cp  =  0.0797  +  (5.556  x  10-5)T,  (17-4  PH  stainless  steel)  (4) 

was  used  with  the  value  of  Tw  at  the  midpoint  of  the  curve  fit.  The  maximum  variation 
of  Cp  over  the  curve  fit  was  less  than  1.2  percent.  The  value  of  density  used  for  17-4 
PH  stainless  steel  was 


5 

7 

9 

13 

17 

25 

41 


w  =  480.0  lbm/ft3 


For  this  test,  the  wall  recovery  temperature  was  selected  as  0.9  To  to  conform  with 
previous  shuttle  heat-transfer  data.  The  heat  balance  becomes: 

dTw/dt 

h(0.9To)  “  b  ^  Q  g-j-Q  -  Tw 

Estimated  uncertainties  for  the  individual  terms  in  Eq.  (5)  were  used  in  the  Taylor 
series  method  of  uncertainty  propagation  to  obtain  the  uncertainty  in  the  heat-transfer 
coefficient  as  given  below: 

h(0.9  To)  Uncertainty,  percent 

10-4  14 

10-3  13 

10-2  11.5 
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These  estimated  uncertainties  are  based  on  a  Tw  of  540“R  and  a  value  of  Tw/Tq 
equal  to  0.818. 


4.0  RESULTS  AND  DISCUSSION 
4.1  EXTERNAL  TANK  ALONE 

During  the  initial  phase  of  the  test,  data  were  obtained  on  the  external  tank  alone. 
These  data  formed  the  baseline  for  comparison  with  the  data  obtained  with  the  integrated 
model  and  provided  an  indication  of  overall  data  quality.  Figure  8  is  a  comparison  of 
0  and  180-deg  ray  heat-transfer-rate  distributions  on  the  external  tank  with  the  smooth 
nose.  Calculated  values  of  laminar  and  turbulent  heating  rates  for  this  configuration  are 
also  presented  in  Fig.  8a  for  zero  angle-of-attack.  The  method  used  to  calculate  these 
heating  rates  is  that  of  Refs.  2,  3,  and  4  with  pressures  calculated  from  modified  Newtonian 
theory.  The  measured  heat-transfer  rates  ahead  of  the  boundary-layer  trip  location,  xj/Lj 
=  0.069,  are  in  good  agreement  with  the  calculated  laminar  rates.  The  trip  effectiveness 
is  evident  with  the  transition  to  the  turbulent  level  completed  at  a  value  of  xt/Lt  of 
0.08.  The  agreement  of  the  measured  and  calculated  turbulent  heat-transfer  rates  is  very 
good  with  the  exception  of  two  regions.  In  the  region  between  Xt/Lt  of  0.2  and  0.35, 
the  heating  rates  indicate  that  the  flow  may  be  relaminarizing.  This  same  trend  was  noted 
consistently  throughout  the  test  data.  An  investigation  of  relaminarization  under 
compressible  flow  conditions  is  reported  by  Nash-Webber  in  Ref.  5.  In  this  work,  a 
boundary-layer  edge  acceleration  parameter  and  a  momentum  thickness  Reynolds  number 
based  on  edge  conditions  are  used  to  correlate  relaminarization. 

Computations  of  these  parameters  for  the  Xx/Lj  location  of  0.25  confirm  the 
probable  occurrence  of  relaminarization  for  the  present  tests.  The  second  region  where 
the  measured  and  calculated  heat-transfer  rates  are  not  in  agreement  is  around  Xj/Lx 
of  0.9.  This  could  be  caused  by  an  upstream  effect  of  the  expansion  at  the  rear  of  the 
tank.  The  tangent  point  for  the  dome  at  the  rear  of  the  tank  is  located  at  Xj/Lj  of 
0.93. 


The  windward  and  leeward  centerline  heat-transfer-rate  distributions  for  the  external 
tank  with  the  smooth  nose  configuration  are  presented  in  Fig.  8b  for  a  10-deg  angle  of 
attack.  The  transition  to  turbulent  heating  is  again  evident  between  Xx/Lx  of  0.06  to 
0.08.  Windward  centerline  heating  rate  measurements  are  in  general  agreement  with 
calculated  turbulent  values.  Note  that  the  lee  side  heating  rates  at  xj/Lj  values  greater 
than  0.6  equal  the  zero  angle-of-attack  case.  This  is  indicative  of  a  vortex  interaction 
phenomenon  as  observed  on  other  configurations  (see  Ref.  6,  for  example). 


10 


AEDC-TR-75-153 


A  comparison  of  axial  heat-transfer-rate  distributions  on  the  external  tank  at  zero 
angle  of  attack  both  with  and  without  the  boundary-layer  trip  is  presented  in  Fig.  9  for 
the  nipple  nose  configuration.  It  is  noted  that  the  distribution  with  the  boundary-layer 
trip  is  essentially  the  same  as  the  clean  distribution.  Upon  closer  examination  it  is  evident 
that  transition  to  the  turbulent  distribution  occurs  between  xt/Lt  of  0.02  and  0.03.  Since 
the  boundary-layer  trip  was  located  at  0.069,  these  data  indicate  that  the  boundary  layer 
is  tripped  by  the  shock  wave/boundary-layer  interaction  generated  by  the  nipple  nose. 

Also  shown  in  Fig.  9  is  a  fairing  of  the  smooth  nose  data  from  Fig.  8a.  It  is  apparent 
that  the  nipple  nose  has  essentially  no  effect  on  the  heating  rates  downstream  of  xt/Lj 
=  0.1. 

Zero  angle-of-attack  axial  heat-transfer-rate  distributions  on  the  external  tank  with 
the  nipple  nose  are  shown  for  Mach  numbers  2.5,  3.5,  4.5,  and  5.5  in  Fig.  10.  The  calculated 
smooth  nose  turbulent  heating  rates  are  included  for  each  Mach  number.  The  agreement 
between  the  measured  and  calculated  values  is,  in  general,  very  good. 

4.2  INTEGRATED  MODEL 

The  axial  heat-transfer-rate  distributions  on  the  leeward  centerline  of  the  external 
tank  for  the  integrated  model  are  presented  for  Mach  numbers  3.5,  4.5,  and  5.5  in  Fig. 
1 1 .  Data  fairings  for  the  external  tank  alone  are  included  to  provide  a  baseline  for 
comparing  the  interference  heating  to  undisturbed  heating  rates.  Peak  heating  resulting 
from  interaction  with  the  orbiter  occurred  at  an  xt/Lt  0.445.  The  ratio  of  the  peak 
heating  to  the  undisturbed  heating  value  is  10  for  Mach  number  3.5  and  5.5  and  12 
for  Mach  number  4.5.  Figure  12  presents  shadowgraph  photographs  of  the  flow  on  the 
external  tank  with  the  nipple  nose  and  on  the  integrated  model,  respectively.  Figure  12b 
shows  the  interaction  of  the  orbiter  bow  shock  with  the  external  tank,  and  the  region 
of  peak  heating  shown  in  Fig.  1 1  is  indicated  on  the  photograph. 

The  heat-transfer-rate  distribution  on  the  side  of  the  external  tank  resulting  from 
the  interaction  of  the  solid  rocket  booster  is  shown  in  Fig.  13.  The  peak  heating  occurred 
at  an  xj/Lj  of  0.35.  The  ratio  of  the  peak  heating  to  the  undisturbed  heating  value 
is  22.  The  bow  shock  from  the  solid  rocket  booster  can  be  seen  in  Fig.  12b. 

5.0  CONCLUDING  REMARKS 
5.1  EXTERNAL  TANK  ALONE 

The  external  tank  was  tested  with  two  nose  tip  configurations  to  obtain  undisturbed 
heat-transfer-rate  distributions.  These  distributions  were  used  to  confirm  the  data  quality 
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and  to  provide  a  comparative  baseline  for  the  heating  rates  from  the  integrated  model. 
The  more  significant  comments  concerning  these  tests  are: 

1.  The  measured  heating  rate  distributions  were,  in  general,  in  good  agreement 
with  calculated  turbulent  heat-transfer-rate  distributions. 

2.  The  nipple  nose  configuration  on  the  external  tank  effectively  tripped  the 
boundary  layer  to  turbulent  conditions  without  the  addition  of  artificial 
roughness. 

3.  A  trend  toward  relaminarization  was  noted  on  the  external  tank  at  xt/Lt 
of  0.25  where  the  nose  ogive  mates  with  the  cylindrical  body. 

5.2  INTEGRATED  MODEL 

The  integrated  model  was  tested  to  investigate  heating  rates  during  the  ascent  phase 
of  the  mission.  Because  of  the  complex  shock  interaction,  no  calculated  heating  rate  values 
were  available  for  comparison.  Instead,  the  peak  heating  values  obtained  on  the  external 
tank  were  compared  with  the  undisturbed  heating  values.  The  peak  heating  rate  on  the 
leeward  centerline  of  the  external  tank  occurred  at  an  Xj/Lt  of  0.44  and  ranged  between 
10  to  12  times  the  undisturbed  value.  No  noticeable  Mach  number  effect  in  the  peak 
heating  was  indicated.  A  peak  heating  ratio  of  22  occurred  on  the  side  of  the  external 
tank,  0T  =  90  deg,  at  an  xt/Lj  of  0.35. 
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X  =  1.125 


a.  Smooth  nose 


b.  Nipple  nose 

Figure  4.  External  tank  nose  tip  configurations. 
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Figure  6.  Photograph  of 
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a.  a  =  0 


b.  a  =  10  deg 

Figure  8.  Heat-transfer-rate  distributions  for  external  tank 
with  smooth  nose. 
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Figure  9.  Influence  of  boundary-layer  trip  on  the  heat-transfer-rate 
distribution  of  the  external  tank  with  nipple  nose. 
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O  External  Tank  with  Integrated  Model 
—  External  Tank  Alone  (see  Fig.  10) 


a.  M^  =  3.5 


b.  M^  =  4.5 

Figure  11.  Heat-transfer-rate  distributions  for  the  externai  tank 
with  interaction  from  the  orbitar. 
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a.  External  tank  with  nipple  nose 
Figure  12.  Model  flow-field  shadowgraphs. 
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b.  Integrated  shuttle  model 
Figure  12.  Concluded. 
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O  External  Tank  with  Integrated  Model 
—  External  Tank  Alone  (see  Fig.  10), 


V^T 


Figure  13.  Heat-transfer-rate  distribution  on  the  side  of  the  external  tank 
with  interaction  from  the  solid  rocket  booster. 
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Table  1.  Test  Matrix 


Model  Configuration 
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0 
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NOMENCLATURE 

Model  skin  thickness,  ft 
Model  wing  span,  in.  (see  Fig.  2) 

Vertical  tail  span  with  origin  at  Z  =  500,  in.  (full  scale)  (see  Fig.  2) 
Specific  heat,  Btu/lbm-°R 

_  q  Btu 

Heat-transfer  coefficient  based  on  To,  ,  ftS.sec-“R 

q  Btu 

Heat-transfer  coefficient  based  on  0.9To,  (q  9To)  -  Tw  •  ft2  -sec-°R 

Reference  heat-transfer  coefficient  based  on  Fay-Ridell  theory  and  a  1-ft 
nose  radius  scaled  to  the  model  scale  (0.0175  ft) 

Axial  length  of  orbiter  model,  in.  (see  Fig.  2) 

Axial  length  of  solid  rocket  booster  model,  in.  (see  Fig.  2) 

Axial  length  of  external  tank  model,  in.  (see  Fig.  2) 

Frce-stream  Mach  number 

Tunnel  stilling  chamber  pressure,  psia 

Heat-transfer  rate,  Btu/ft^-sec 

Free-stream  unit  Reynolds  number,  ft'^ 

Tunnel  stilling  chamber  temperature,  °R 

Model  wall  temperature,  °R 

Time,  sec 

Model  skin  density,  Ibm/ft^ 

Longitudinal  coordinate  with  origin  at  the  orbiter  model  nose.  in.  (see  Fig. 

2) 

Longitudinal  coordinate  with  origin  at  the  solid  rocket  booster  model  nose, 
in.  (see  Fig.  2) 
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Xt  Longitudinal  coordinate  with  origin  at  the  external  tank  model  nose,  in. 

(see  Fig.  2) 

Y  Lateral  coordinate,  in.  (see  Fig.  2) 

Z  Vertical  coordinate,  in,  (see  Fig.  2) 

a  Angle  of  attack,  deg 

0  Angle  of  sideslip,  equal  to  negative  yaw  angle,  deg  (see  Fig.  2) 

0o  Angular  measurement  on  orbiter  model,  deg  (see  Fig.  2) 

0s  Angular  measurement  on  solid  rocket  booster,  deg  (see  Fig.  2) 

0T  Angular  measurement  on  external  tank,  deg  (see  Fig.  2) 

SUBSCRIPTS 

i  Initial  conditions 
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